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In this work, we investigate the 4d path integral for Euclidean quantum gravity on a hypercubic 
lattice, as given by the Spin Foam model by Engle, Pereira, Rovelli, Livine, Freidel and Krasnov 
(EPRL-FK). To tackle the problem, we restrict to a set of quantum geometries that reflects the 
large amount of lattice symmetries. In particular, the sum over intertwiners is restricted to quantum 
cuboids, i.e. coherent intertwiners which describe a cuboidal geometry in the large-j limit. 

Using asymptotic expressions for the vertex amplitude, we find several interesting properties of 
the state sum. First of all, the value of coupling constants in the amplitude functions determines 
whether geometric or non-geometric configurations dominate the path integral. Secondly, there 
is a critical value of the coupling constant a, which separates two phases. In both phases, the 
diffeomorphism symmetry appears to be broken. In one, the dominant contribution comes from 
highly irregular, in the other from highly regular conhgurations, both describing flat Euclidean 
space with small quantum fluctuations around them, viewed in different coordinate systems. On 
the critical point diffeomorphism symmetry is nearly restored, however. 

Thirdly, we use the state sum to compute the physical norm of kinematical states, i.e. their norm 
in the physical Hilbert space. We find that states which describe boundary geometry with high 
torsion have exponentially suppressed physical norm. We argue that this allows one to exclude 
them from the state sum in calculations. 


I. MOTIVATION 

The spin foam approach has been developed to give 
a rigorous meaning to the path integral for quantum 
gravity (see [l[ for a review). Its central idea rests on 
the observation that the first order formalism of GR can 
be rewritten as a certain constrained topological theory 
Hi, dubbed “BE theory”. A quantization choice for 
these constraints is what specifies the spin foam model, 
and in recent years there have been several proposals 
A popular choice has emerged in the so-called 
EPRL-FK model @,10, which possesses quite useful 
properties. In particular, the resulting amplitude has an 
asymptotic expression for larg e q uantum numbers which 
reproduces the Regge action Furthermore, the 

resulting path integral for the EPRL-FK model natu¬ 
rally has boundary states which resemble the spin net¬ 
work states from canonical loop quantum gravity [iMl, 
which is why it has been coined “covariant loop quantum 
gravity” ([l6|, see also Haiii). 

There are several open questions, however. While there 
are numerous results available which elucidate the prop¬ 
erty of a single vertex amplitude, very little is known 
about the behaviour of the whole path integral Q In par¬ 
ticular, it is the realm of many building blocks which 
is of utmost importance if one wants to understand the 
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^ There are some results on the asymptotic expression for more 
than one vertex [T^.[^. as well as self-energy calculations [^. 


continuum limit of the theory. This is closely connected 
to the question of renormalizability of the model, as well 
as the challen ge o f restoring the broken diffeomorphism 
invariance |22l-[^ , both of which are unsolved questions 
up to this point. 

The non-perturbative and background-independent 
nature of loop quantum gravity suggests that numerical 
methods from lattice gauge theory could be very useful 
tools for investigations, as well as extracting physical pre¬ 
dictions from the theory. This assumption is supported 
by the fact that spin foam models are by construction 
generalized lattice gauge theories Hi 

There has been considerable progress in recent years on 
understanding the continuum limit of so-called spin net 
models, which are analogues of spin foam models [28|. 
Using numerical algorithms in tensor network renormal¬ 
ization HH, many results about the phase structure 
of these models, as well as their continuum limit proper¬ 
ties, can be derived, in p articular in the realm of finite 
and quantum groups [^, . The methods can also ap¬ 

plied to lattice gauge theories [33| . 

This, together with the fact that many results in lat¬ 
tice gauge theory are accessed numerically, in particular 
high-precision predictions of physical quantities, suggests 
that the numerical investigation of the full spin foam path 


^ A complementary but closely connected road towards the un¬ 
derstanding of the continuum limit comes from the group field 
theory / tensor field theory approach imi, which is using 
much more the particle physics understanding of the spin foam 
amplitudes. 
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integral will be a crucial step on the road to a predictive 
quantum gravity theory, as well as in understandin g it s 
continuum limit and renormalization group flow [3, [35| . 
It is the goal of this article to make a step into this di¬ 
rection. 


The symmetry-restricted state sum 

In many calculations performed - both analytically 
and numerically - in lattice gauge theory, the art is to 
approximate the path integral in the right way. The goal 
is to make it simple enough to handle, while still keeping 
an expression from which physical information can be ex¬ 
tracted. One way to perform this is to not sum over all 
histories of the fields in the path integral, but only over 
those which are relevant for the physical process one is 
interested in. That way, one can still hope to gain realis¬ 
tic approximations for the values of certain observables. 
Of course, the approximation in question therefore has 
to be tailored to the situation one wants to investigate. 

We intend to copy this strategy for the spin foam ap¬ 
proach in this article. We believe that this could, if done 
in the right way, provide very useful approximations for 
the spin foam state sum. Hence this might ultimately 
provide a pathway to statements not only about the con¬ 
tinuum limit of spin foams, but about expectation values 
of physical observables, as well as the renormalization 
group flow of the model. In all of these, the spin foam 
state sum has to be performed, which in full generality 
can most likely not be achieved. 

As a first step, in our article we will look at a drastic 
approximation of the whole path integral. These will ren¬ 
der the physical predictability of the model questionable 
in certain regions. However, the model will still be com¬ 
plicated enough so that certain concepts can be tested, 
and specific open questions can be investigated. 

In particuiar, we wiii work on a regular hypercubic lat¬ 
tice in 4d. Furthermore, on this iattice, we wiii not con¬ 
sider aii possibie states, but only those which conform to 
the lattice symmetry. This is a condition on the inter¬ 
twiners, which we pick to correspond to cuboids. 0 

A cuboid is completely determined by its three edge 
lengths, or equivalently by its three areas (see figure [1]). 
Ail internal angles are and the condition of regular 
cuboids on all dual edges of the lattice will result in a 
high degree of symmetries on the labels: The area (and 
hence the spin) on each two parallel squares of the lattice 
which are translations perpendicular to the squares, have 
to be equal. 

A few comments about this simplification are in order: 


^ This choice has also been considered in [3^ in the canonical con¬ 
text, and in [ 3 ^ in the cosmological setting. The hypercubic 
lattice has also been investigated in the loop quantum gravity 
path integral framework in l40l| . 
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FIG. 1. A cuboid in R^. The face normals rii are orthogonal, 
and normalized by the face areas. 


• The main reason for this symmetry is that it sim¬ 
plifies the problem dramatically, making the sum¬ 
mation over states much more manageable. Also, 
the lattice is directly tailored for making contact 
with usual forms of lattice gauge theory, which we 
will explore by coupling matter degrees of freedom 
in a companion paper, in the spirit of 0. 

• As one can readily see, the high degree of symmetry 

will make all quantum geometries flat. The analysis 
carried out in this article is therefore definitely not 
suited for describing local curvature. Still, there 
remain enough states to describe local degrees of 
freedom. These will be non-geometric^ as a result of 
the twisted geometries appearing in the spin foam 
formalism [42|, • The role and behaviour of this 

non-geometricity will form a major part in this ar¬ 
ticle. 

• The set of considered states is still large enough so 
that an infinite-dimensional (Abelian) subgroup of 
the diffeomorphisms acts on it. This makes this 
model an ideal test-bed for the question of how to 
treat these diffeomorphisms on the lattice, and how 
to treat them in the state-sum. 

• Restricting the state sum to only a subclass of 
states is not only a tool to simplify the state sum. It 
results in an approximation which can be valid for 
describing situations in which the removed states 
are not suspected to have a great influence on the 
physics. For instance, one could disregard states 
with high curvature if one is only interested in the 
low-curvature regime of the path integral. While 
considering only quantum cuboids is certainly a 
very severe restriction, one could try to access, say, 
the propagation of gravitational waves on a lattice, 
by adding very specific states describing low local 
curvature. A controlled enlargement of the states 
beyond the set of quantum cuboids could therefore 
provide an excellent tool for approximations of the 
spin foam path integral. 

The plan of our article is as follows: 
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We will give a brief review of the EPRL-FK spin foam 
model in section mi Before performing computations in 
the quantum regime, we will look at the semiclassical 
regime of the path integral. The corresponding geomet¬ 
rical intuition, which we will develop in section mil will 
be very useful for the rest of the article. 

The actual construction of the quantum cuboid inter¬ 
twiner will be carried out in section lEl Afterwards, the 
full vertex amplitude will be considered in section m 
where we in particular describe its asymptotic expres¬ 
sion for large spins (i.e. areas). Finally, the results of 
our numerical investigation of the quantum path integral 
for more than one building block, in particular in terms 
of the geometric interpretation developed earlier, will be 
given in section lYD 

We will see that, despite the severe simplification of the 
system by restricting to quantum cuboids, the analysis 
will reveal a surprising amount of insight into the path 
integral. 

II. INTRODUCTION 

The spin foam state sum we will employ in this article 
will be the Euclidean EPRL-FK model with Barbero- 
Immirzi parameter 7 < 1. The main reason for this is 
that this model is, at the current state, one of the best 
understood ones, with a lot of literature already existing 
on the subject. 

The EPRL-FK model was originally defined on a 2- 
complex dual to a 4d triangulation, but there is a clear 
generalization to arbitrary 2 -complexes [i^, which we 
will use in this article. A 2 -complex T is determined by 
its vertices its edges e connecting two vertices, and 
faces / which are bounded by the edges. While this data 
is parley combinatorial, we will think of T as being em¬ 
bedded in a manifold. For practical purposes one needs 
to choose fiducial orientations of edges and faces, where 
however the result does not depend on that choice. 

The spin foam model can either be written in an oper¬ 
ator formulation , or in a dual holonomy formulation 
[^ . In the former one, the path integral is formulated 
as a sum over states. A state in this context is given 
by a collection of spins, i.e. irreducible representations 
jf ^ of SU{2) to the faces, as well as a collection of 
intertwiners Lq on edges. Here an intertwiner on the edge 
e is an invariant tensor in the space Vj^ 0 ... 0 Vj ^, 
where ji,... jn are the spin on the faces meeting at e, and 
the representation spaces are either the usual one or 
the dual representation space for that spin, depending 
on whether the relative orientations of e and / agree or 
disagree. 

The actual sum is given by 

j/,^e f e V 

where Af, Ae and Ay are the faee-, edge- and vertex- 
amplitude functions, depending on the state. The sum 


has to be carried out over all spins, and over an orthonor¬ 
mal basis in the intertwiner space at each edge. 0 

The allowed spins jf in the EPRL-FK model are such 
that := both are also spins, i.e. half-integers 0 

Popular choice^ for the face amplitudes are either 
Af = 2jf + 1 or = ( 2 jf^ + + 1). In what 

follows, we keep an open mind about the precise form 
of the face amplitude, and allow for some freedom. We 
in particular introduce an additional parameter a, and 
write 

((27 +1) (27+1))“, (2) 

and investigate the influence different choices for a have 
on the state sum. 

The edge amplitudes Ae are usually taken to be equal 
to 1 in our formulation, but there are certain ways of 
formulating the vertex amplitude, in which one puts the 
normalization of the intertwiners, or left over sign factors 
here. 

The main ingredient to the state sum o is the vertex 
amplitude, the construction of which we briefly describe 
in what follows. For this we need to describe the boost¬ 
ing map 4>, which maps the 5'I/(2)-intertwiners te to an 
SU (2) X SU (2)-intertwiner by 

$ : Invs[/(2) Vj. — Invspi„(4) V-+ 0 Vj- . 

i i 

The map 4> is given by 4 > = 7 ^ o where V is the 
projector on SU{2) x S'I/(2)-intertwiners, and the map 
f3 : Vj ^ Vj+ ^ Vj- is given, in the case that 7 < 1 
which we use throughout this article, by the isometric 
embedding of Vj into Vj+_^j_, the highest weight space 
in the Clebsh-Gordan decomposition of Vj+ G Vj- in the 
SU(2) representation category. 

The vertex amplitude at the vertex v is the contraction 

Av = tr(^(^$(;,e)) (3) 

eDv 


There is a subtlety at this point: if one sees the state sum as 
coming from a restriction of the 4d BF amplitude on which the 
simplicity constraints have been imposed, one should actually 
not sum over an orthonormal basis of S't/( 2 )-intertwiners, but 
rather over a set Le such that 4>(6e) form an orthonormal subset 
of the SU (2) X SU (2)-intertwiners, where is the boosting map. 
These two choices are not the same, since it can be shown that 4> 
is not an isometry [^. Since in this article, we will only have one 
specific intertwiner for each edge, the issue does in fact not arise. 
It will, however, emerge as soon as one allows for intertwiner 
fiuctuating in the path integral. 

^ This means that, depending on 7 , the set of allowed spins can 
actually be very small, or even just be the trivial representation, 
if 7 is irrational. This feature is very particular to the Euclidean 
version of the model, and disappears in the Lorentzian setting, 
where all real 7 are allowed. 

® This choice has infiuence on the convergence of the state sum, 
see e.g. [21I. 
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where one takes either 4>(^e) or the dual (4>^e)^ depend¬ 
ing on whether the edge e is incoming or outgoing of 
the vertex v. The tr-contraction is to be understood in 
the sense that the tensor product in ([3]), for each face / 
touching there will be two 5'pm(4) indices in the rep¬ 
resentation which can be contracted since they 

appear in opposite position. 


A. Coherent intertwiners 

A particular advancement in the spin foam model ap¬ 
proach was the observation that the intertwiners te can 
be given a geometric interpretation in terms of polyhedra 
in |4§-[^. For a spin j, denote by \j j) the highest 
weight vector with respect to rs = |cr 3 , corresponding to 
the unit vector in the 3- or 2 ;-direction. For any other unit 
vector n G 5'^, choose g G SU{2) such that g > = n, 

then the coherent vector 


\jn) = g\>\jj). 


(4) 


Given a collection of spins ji,... jn and vectors ni,... 
which close, i.e. = 0, one can define the coherent 

polyhedron 


I'') = / dg g>(^\jifii). 

Jsu( 2 ) 7 


ISU{2) 


(5) 


The geometric interpretation is that of a polyhedron, 
with face areas jf and face normals The closure 
condition ensures that such a polyhedron exists, and is 
unique up to translation. The formula ^ also exists 
for non-closing normals, which do not have a polyhedral 
interpretation. However, it can be shown that the state 
sum model m can be written solely as an inte gral over 
coherent polyhedra, i.e. having closing normals |5l| . 

This concludes the recap of the general setup. In 
the following we will describe the specific simplifications 
within this setup that we will use in order to obtain phys¬ 
ical intuitions for the full path integral. 


later be encoded in the state, by specification of spins jf 
and intertwiners Lq. 

In what follows, we consider the geometric interpreta¬ 
tion of the semiclassical regime, not only of one inter¬ 
twiner, but of the whole lattice. We are in particular 
interested in the large j-regime of the quantum cuboids. 
In this limit, these become classical cuboids as in figure [TJ 
which are completely specified by their three areas. The 
fact that opposite areas are identical lead to a transla¬ 
tion symmetry in the areas of the lattice. Therefore, a 
semiclassical configuration is given by an assignment of 
areas a = jTp to the squares of the hypercubic lattice, 
where two areas agree whenever they are parallel and one 
is reached from the other by a translation perpendicular 
to them. 

Denote the four directions in the lattice by 
and numerate the hypercubes by n G The areas 
then satisfy 


n+pcp+qecr _ 


(6) 


with all different directions, ep^Ca- unit vectors 

in p, and a direction, and p^q ^7L. 

Dual to one vertex n in the lattice is then a semiclassi¬ 
cal quantum hypercuboid^ which is determined by the val¬ 
ues of six of its areas ..., . Since a classical 

hyper cuboid in R^ is determined by its four edge lengths, 
this means there are two excess degrees of freedom. It 
is straightforward to see that these cannot correspond 
to geometric configurations: Assume that six areas are 
given, these determine uniquely the three areas for each 
individual 3d cuboid in the hypercuboid’s boundary. In 
turn, for each such 3d cuboid, the three areas determine 
three edge lengths. In order for there to be a geomet¬ 
ric hypercuboid, not only do the areas of two cuboids 
meeting in a rectangle have to match, but both individ¬ 
ual edge lengths of that rectangle have to agree. It is 
not difficult to derive that the conditions the areas have 
to satisfy, so that both edge lengths of each rectangle, 
as seen from neighboring cuboids, agree. This condition 
can be written as 


III. SEMICLASSICAL CONSIDERATIONS 


^xy^zt — ^xz^yt — ^xt^yz- 


(7) 


In this article, we will consider the case where the 
space-time manifold J\4 ^ x [0,1] is the product of the 
3-torus and a closed interval. Hence space is compact- 
ified toroidally. We cover M. by dd hyper cubes, which 
form a regular hypercubic lattice H. The 2-complex F 
we use is the 2-skeleton of the dual complex. So there 
is a vertex for each hyper cube, and two vertices are con¬ 
nected by an edge whenever two hypercubes intersect in 
a 3d cube. The faces of F are dual to squares in H, on 
which four hypercubes meet. 

Note that this hypercubic lattice does not carry any 
geometric information at this point. The geometry will 


It seems worthwhile noting that this problem - of match¬ 
ing areas but non-matching shape - is well-known, and 
usually referred to as twisted geometries in the literature 


The 4-volume of a geometric hypercube is given by 
the product of areas of two perpendicular faces, so the 
geometricity constraints © are the condition that the 
three choices of opposite pairs among the six faces deliver 
the same result. For a non-geometric configuration, we 
define the 4-volume of a hypercube as 

1 

V •— {f^xy^xz * ‘ ‘ ? 


(8) 


which reduces to the original definition at geometricity. 
If one were to define the four diameters to be 
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1^4 


(ayzaytatz)'^ 


if^xy^xz^xt) ^ 
{ayzaytatz)^ 


(9) 


with similar definitions for the other three directions in 
the lattice, one would have that V 4 = dxdydzdt, even for 
non-geometric configurations. We also define the non- 
geometricity 


r := 


( ^xy^zt ^xz^yt 
(^xz^^yt ^xt^yz 
(^xy^zt ^xt^yz 


( 10 ) 


as a measure of the deviation from the constraints O- 


A. Lattice deformations 
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FIG. 2. Action of a diffeomorphism in t-direction on the areas 


The fact that the hypercubic lattice in the manifold is 
not prescribed with an ab initio geometry, but obtains ge¬ 
ometric interpretation via the state, has important con¬ 
sequences. In particular, there are some configurations 
which arise as lattice deformations of one another. These 
can be seen as the result of either an active, or a passive 
diffeomorphism. 

Imagine, for definiteness, the manifold to be equipped 
with a standard, fiat Euclidean metric 77 , and the lat¬ 
tice T inserted being completely regular, i.e. all lengths 
(and hence all areas) are equal. Now consider a diffeo¬ 
morphism (j) which deforms the metric in a thin 3d strip, 
say, the slightly thickened hyperplane around t = 0 . As¬ 
sume that (j) is such that it stretches space below, and 
contracts the space above the hyperplane. As a result, 
with respect to the new metric g = ((/)“^)* 77 , the areas 
of the “time-like” squares above the hyperplane will be 
smaller, and the ones below it will be larger (see figure 

Ej). 

So this diffeomorphism of the manifold M. induces a 
change in the areas of the lattice, since it changes the 
metric, but leaves the lattice where it is (if it were to 
transform the lattice in a similar way, the areas would 
remain unchanged trivially). This can be regarded as 
an “active” diffeomorphism - the ’’passive” version can 
be seen as keeping the metric fixed, but changing the 
embedding of the lattice T ^ 0 (r), by moving the t = 0 - 
hyperplane into positive t-direction. The result on the 
values of the areas in the lattice will be the same, so 
by just looking at the lattice and not at the manifold, 
the difference between “active” and “passive” diffeomor¬ 
phism vanishes. 

The continuum diffeomorphisms prescribed in this way 
are generated by vectorfields which only act in t direc¬ 
tion on the manifold, i.e. X(x, 2 ), t) = T{t)dt^ with T a 
smooth function. Together with the deformations in the 
other three major directions, the resulting collection of 
diffeomorphisms forms an Abelian subgroup of the full 
group of Diff(AI). 


The variable transformation we have just described can 
be extended to the possibly non-geometric configurations 
that arise in the large-j-regime of the EPRL-FK model. 
Just as in the example, there is one generator of the group 
for each 3d hyper surface in the lattice, “moving” that hy¬ 
persurface in a direction orthogonal to it, while keeping 
the rest of the lattice fixed. As one can see, this group 
is finite-dimensional for compact manifolds, but the di¬ 
mension grows as the number of hyper cubes increases, 
becoming infinite in the limit of infinitely large lattices. 

To describe the action of these “lattice deformations” 
on the areas consider, for definiteness, a 3d hyper sur¬ 
face in the x^z-directions. Let n be a hypercube whose 
future xyz cuboid is part of the chosen hypersurface, then 
the transformation is defined to be 

+ i = x,y,z (11) 

for all such n, while remains unchanged. The time- 
like areas of the following hypersurface gets changed as 
well, by 

^afA‘(l-4"+"‘) i = x,y,z. ( 12 ) 

The parameters A’^, quantify the stretching 

and compressing of the respective hyper surfaces. It 
is easy to see that, in order not to violate A" 
has to be the same among the whole hyper surface, 
i.e. = A^ =: Ai for all m^. A simi¬ 

lar statement holds for A 2 := A’^+^E 

There is an interesting observation to be made at this 
point: In general a diffeomorphism should not change the 
total 4-volume of space-time. Depending on the areas 
there might be no choice for A 1^2 so that the 4-volume 
is preserved, however. So, while the action of the diffeo 
group can be defined on all configurations, there are only 
some on which it can be defined in a 4-volume-preserving 
way. We call these configurations semi-geometric. A suf¬ 
ficient condition for semi-geometricity is that the sum of 
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the 4-volumes of a pair of hypercubes n, n + stays con¬ 
stant under a diffeomorphism, so 4-volume is not only 
preserved in general, but also locally. It is not difficult 
to show that this is equivalent to the following condition: 
For every 2 x 2-arrangement of squares in the lattice, the 
two respective products of diagonally opposite areas in 
that arrangement are equal (see figure [3j). It should be 
clear that the diffeomorphisms, as we have defined them 
here, preserve the notion of semi-geometricity. 


ai 

a2 

as 



FIG. 3. The areas of flat squares in Euclidean space arranged 
in a tile satisfy aia 4 = a 2 a 3 

It should also be noted that semi-geometricity is 
weaker than geometricity ©: a totally translation- 
symmetric configuration where all hypercubes are con¬ 
gruent, is certainly semi-geometric. However, the six ar¬ 
eas of the hypercube can be chosen in a non-geometric 
way. Note that geometricity is purely local, in terms of 
one hypercube, while semi-geometricity is a condition on 
relations of areas between neighboring hyper cubes. 


eight cuboids (see figure |4)). The vertex amplitude is de¬ 
fined with the help of a boundary spin network F, which 
in our case is of the special form depicted in figure [5] . 
The amplitude is given in terms of spins ji associated to 
links I of F, as well as coherent intertwiners La © asso¬ 
ciated to the nodes a of F. The intertwiner can be given 
in terms of normal vectors which is thought of be¬ 
ing the normal vector of the link I = (a 6 ), sitting at the 
polyhedron at a. 



FIG. 4. A four-dimensional hypercuboid represented by its 
boundary, consisting of eight cuboids. They are labelled by 
a = 0 ... 7, as depicted. 


IV. QUANTUM CUBOIDS 


We now turn to the quantum theory. 

In the 2-complex we consider, every edge has six faces 
attached to it, corresponding to the six faces of the cubes. 
So any intertwiner in the state-sum will be six-valent, 
and therefore can be described by a coherent polyhedron 
with six faces. In our setup, we restrict the state-sum 
to coherent cuboids, or quantum cuboids. A cuboid is 
characterized by areas on opposite sides of the cuboid 
being equal, and the respective normals being negatives 
of one another (see figure [ 1 ]). 

The state Gi,i 2 ,i 3 given by 


= / dg g>(^\ji,ei)\ji,-ei). (13) 
JsUi2) 

Here ei = exp(—i 7 rcr 2 / 4 ) > 63 , 62 = exp(i 7 rcri/ 4 ) > 63 , and 
63 are taken to be unit vectors in R^. It is worth noting 
that the intertwiner space is nonempty for every choice 
of spins ji,^ 2 , 43 , and the state (p!3|) always exists and is 
non-vanishing. 



FIG. 5. The spin network for a hypercuboid. Each node is six 
valent, and spins of edges emanating opposite of each other 
from a node coincide. Therefore, only six different spins exist. 
In the semiclassical limit these become the six areas of 
the hypercuboid. 


V. THE AMPLITUDE 

The vertices in the four-dimensional hypercubic lattice 
are all eight-valent, since a 4d hypercube is bounded by 


The vertex amplitude for a Barbero-Immirzi parameter 
7 < 1 factorizes as Ay = Ay A ~, with 

At = [ dga (14) 

J SU{2)^ 
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with the complex action 

= t^^2j;ln(-nab|5(^^5(6|nba) (15) 

I 

=■■ -^S[ge] 

where, in formula dm, a is the source node of the link 
/, while b is its target node, and all states are in the 
fundamental representation. 


A. Large j asymptotics 

The amplitudes possess an asymptotic expression 
for large ji , which has been investigated in [ 10 , • To 

compute this expression, we follow the analysis in la, 
and refer to that article for calculational details. 

Firstly, we note that out of the N = S group integra¬ 
tions in (HI, one is obsolete due to the invariances of 
the Haar measure dg. The remaining integral possesses 
isolated critical points, and can therefore be evaluated by 
an extended stationary phase approximation. The action 
m is readily seen to be invariant under any ga —ga^ 
so there is a 2 ^-fold symmetry. 

Modulo this symmetry, there are two distinct station¬ 
ary and critical points, satisfying the equations 

nab ’= Qa t> flab = -Qb > ^6a (16) 

for all links {ah) (the closure condition is automatically 
satisfied by the choice of quantum cuboids, for any 
selection of spin ji). Using the convention in figure IH 
having fixed go = the two solutions Si and S 2 are 
shown in table [H 



El 

E2 

91 

exp (if < 7 i) 

exp (-if< 7 i) 

92 

exp (if ( 72 ) 

exp (-if <72) 

93 

exp (if (73) 

exp (-if <73) 

9^ 

exp (-if (73) 

exp (if <73) 

95 

exp (-if (72) 

exp (if <72) 

96 

exp (-if ( 7 i) 

exp (if< 7 i) 

97 

exp (if ( 7 i) 

exp (if < 7 i) 


TABLE I. The two solutions for ,... gr of (cai not connected 
by symmetry, corresponding to the stationary and critical 
points of the action m- Note that ^ 1 ,..., all correspond 
to rotations by 7r/2, while gr corresponds to rotations by tt. 

It can be easily checked that these are, in fact solutions 
to the extended stationary phase approximation equa¬ 
tions (|T6l). Also, one can readily see that these are the 
only solutions: the six equations for the edges ( 01 ) up 


to (06) force each of the respective six group elements 
gi {i = 1 ,... 6 ) to lie in some 1 -parameter subgroup, 
corresponding to a rotation around the respective nor¬ 
mal fioi = fioi of cuboid 0. These angles are fixed by 
the equations m for a, 6 G 1,... 6 , and there are two 
overall solutions (modulo 2^ signs). Remaining are six 
vector equations for one SU{2) element ^ 7 , so the sys¬ 
tem appears drastically over deter mined. However, one 
finds that the solution gj = ±exp(i^cri) solves all of 
these simultaneously, and these are the only ones which 
do so. This leaves us with two solutions (up to 2^ signs), 
depicted in table [H 

Another comment about these solutions are in or¬ 
der: As in the case for the 4-simplex, there are “non¬ 
geometric” boundary data, which cannot be used to de¬ 
fine a continuous 3-metric on all of the boundary of the 
hypercuboid (in the 4-simplex case, these are called “non- 
Regge-like”). Unlike in the 4-simplex case, however, 
not all of these non-geometric boundary data are sup¬ 
pressed by the stationary phase approximation method. 
Rather, there are non-degenerate, but non-Regge-like 
states which appear in the large j-limit. These are de¬ 
scribed in section Hm as states violating the geometricity 
condition ©• This issue can be traced back to the dif¬ 
ferences in geometry of a 4-simplex and a 4-hypercuboid, 
in particular the fact that specifying all areas (satisfying 
certain inequialities) in the former determines a geomet¬ 
ric state, while it does not in the latter. This point cer¬ 
tainly deserved more attention, and we will come back to 
it at another poinfQ. 

Each critical stationary point gc contributes one term 
of the form 


I = 


(27r 


i21 


det{-H{ga)) 


oS{gc) 


(17) 


to the approximation formula, where H{gc) is the Hessian 
matrix of S', evaluated on the points g = gc. 

One can readily see that both critical stationary points 
listed in table [I] lead to S[gc] = 0 . This is not surpris¬ 
ing, since all dihedral angles equal and at each square 
in the lattice precisely four hypercuboids meet, leading 
to vanishing deficit angle. In other words, the choice of 
quantum cuboids as the only allowed intertwiners pre¬ 
vents local curvature (intrinsic and extrinsic) to be ex¬ 
cited. Hence the Einstein-Hilbert action part vanishes, 
as well as the boundary terms. 

The stationary phase approximation for the amplitude 
is therefore completely determined by the determinant of 
the Hessian matrix H of second derivatives of S. This 
determines the (first order) quantum corrections to the 
classical solution in the path integral. It therefore con¬ 
tains nontrivial quantum information about the path in¬ 
tegral measure, beyond the Einstein-Hilbert action. 


^ We thank the anonymous referee who pointed out the importance 
of this issue. 
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Introducing, around each critical point coordinates 
{X[} via Qa = ^a,c exp(icr/X^/ 2 ), one obtains for the 
entries of the Hessian matrix 

= - -'^ibnib) (18) 

{ab)Da 

= -i^iJKfiab-nibnib) ( 19 ) 

where the vectors n are given by m- Note that the sec¬ 
ond expression m is symmetric, because fiab = —^ba- 
The Hessian matrix is 21 x 21, and its determinant can 
be computed by a computational algebra program. The 
result is a quite long expression, which can nevertheless 
be handled analytically and numerically. As it turns out, 
the two determinants for Hi and H 2 are complex conju¬ 
gate of each other, so that the whole amplitude is real. 


d^s 

dxidxJ 

d^s 

dXidXJ 


z 


1-7" 

4 


aF-^E+^V 

EHif 


jf f 



In the last line, we have defined the dressed vertex am¬ 
plitude in which the face- and edge amplitudes have 
been absorbed in such a way that boundary amplitudes 
are taken care of correctly. Note that the only way 
the asymptotic state sum (j22j) depends on the Barbero- 
Immirzi paramter 7 is via the prefactor, which does not 
influence the physics. This is a direct consequence of the 
fact that no curvature degrees of freedom are excited in 
our simplified state sum. 


B. Properties of the amplitude 

For the following article, we will use the following no¬ 
tation: the (dressed) vertex amplitude Ay depends on six 
spins jf'i, • • • 5^65 so we write 

^v{j) = j2, is, je), (23) 

using the notation in figure[5]. In particular, ji, 2,3 denote 
the spins on space-like, ^ 4 , 5,6 the spins on time-like faces. 
On some occasions in what follows, we will treat space- 
and time-like faces separately, so we also introduce the 
short-hand 

‘^'^(j:^) • (jl 5 J2 5 J3 5 5 ^2 5 ^2 ) : (24) 

for the amplitude of the equilateral hypercuboid. So 
whenever the amplitude has two arguments, they refer 
to space- and time-like spins, respectively. 


It is worth noting at this point that, because S' = 0 on 
critical stationary points for the quantum hypercuboid, 
the only dependence on the Barbero-Immirzi parameter 
7 is via the Hessian matrix determinant det H. Since 
detH'(Aj) = det H'(j), i.e. H is homogenous in the 
spins, one can readily see that the amplitudes A^ satisfy, 
in the large j limit, 

(^) ' = (^) ' 


Similarly, in the large j-limit, the norm squared of the 
quantum cuboid states (USD is given by 


II 


2 


8(1-7^)-^ 

(ii +i2)(j2 + j3)(ii +h) 


( 21 ) 


With this, we get for the state sum, in the large-j limit 
on a regular hypercubic lattice, that 


n(ji +•^2)(.?2+ h){h +is) n 


( 22 ) 


In the hypercubic lattice the relation 4^ = = "f 

between the numbers of vertices, faces and edges, holds. 
A term Ay in the state sum therefore scales with the 
spins as 


A(AJ) ~ (25) 

This means that for a critical value of a^o = |, the ampli- 
tude becomes scale-invariant. For larger a, the state sum 
is most definitely divergent, while there can be found a 
value of a small enough so that the state sum converges. 

The situation is different if we do not consider the state 
sum for a closed manifold, but one with boundary. Con¬ 
sider a finite regular hypercubic lattice decomposition of 
M = X [0,1]. The spatial boundary manifolds are 
three-tori, and fixing boundary spin networks is equiva¬ 
lent to fixing all spatial spins throughout the lattice, due 
to the symmetry (| 6 ]). The scaling of an amplitude with 
respect to the time-like spins goes like 

XilXk) - A®“-®. (26) 


This means that, when computing the physical inner 
product with the help of the state sum (| 2 ^ on a hypercu¬ 
bic lattice, one can expect the sum to converge absolutely 
for Of < |, when the lattice is large enough. 
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VI. RESULTS 

In the following, we will investigate the asymptotic am¬ 
plitude for the hyper cuboid. This will give an insight into 
the quantitative and qualitative behaviour of the spin 
foam state sum. 


A. Non-geometricity in the path integral 


One interesting question is how much configurations 
with high non-geometricity ^ m contribute to the path 
integral (j22j). To investigate this, we consider a regular 
geometric hypercube, i.e. with all spins equal. 

Then we consider the value of the amplitude Ay , depend¬ 
ing on the two independent non-geometric directions of 
configurations ^ around the regular configuration. We 
vary the three time-like areas a^t, i = x, y, z, by 




^6 


1 

71 

1 

Vg 


(^ScLxt H“ • 


(27) 


Note that this variation is isochoric, i.e. satisfies SV = 0. 
The value of the amplitude is shown in figure [6l 



FIG. 6. Value of the dressed vertex amplitude depending on 
the non-geometricity of the labels. 

The amplitude shows a clear maximum at the geomet¬ 
ric configuration. It is noteworthy that this is not a local 
maximum in all of state space, though. The reason is that 
there are four more directions in which the amplitude can 
be varied, which can change the shape of the hyper cube. 
These directions preserve geometricity. Of those will be 
one direction corresponding to simply scaling the hyper¬ 
cube. Since the amplitude scales as ([26]) with respect 
to this direction, the statement can therefore be phrased 


as follows: Among all configurations with a fixed four- 
volume, the ones contributing most to the path integral 
will be the ones with vanishing non-geometricity. 

The analysis suggests that fluctuations in the non- 
geometricity ^ around ^ = 0 do contribute to the path 
integral, proportional to the width of the Gaussian in fig¬ 
ure El This width m| can be interpreted as a mass term 
for since it measures the ease with which ^ can be ex¬ 
cited locally. The precise value of this mass depends on 
the parameter a in the path integral. The connection 
between the two can readily be computed from the exact 
form of the large-j asymptotic formula, to be 

mlia) = 2a - (28) 

SO for a < III « 0.4854, the mass m| < 0 will be 
negative. This will result in the main contribution to 
the path integral coming from highly non-geometric 
configurations. One could argue that this restricts the 
desired value of a on physical grounds, since the path 
integral should deliver geometric states in the classical 
limit. 

The presence of a mass term suggests that there could 
also be a kinematical term for turning the non- 
geometricity into a propagating degree of freedom. To 
investigate this further, one would not only use more than 
one building block, one would in particular have to con¬ 
sider more intertwiners than just the quantum cuboid, to 
investigate realistic propagation of local excitations of 

In particular, there is the possibility that there ex¬ 
ists a regime in the full path integral in which the non- 
geometricity ^ becomes effectively described by a fluctu¬ 
ating quantum field on some background geometry. This 
point warrants further investigation. 


B. Vertex-displacement symmetry and the 
amplitude 

We now come to an investigation of the amplitude for 
more than one building block. To this end, we first con¬ 
sider two hyper cubes, which are glued together along a 
common (space-like) cuboid. One hypercuboid can there¬ 
fore be regarded as being “in the future” of the other. We 
assume that the main contribution comes from the geo¬ 
metric configurations, i.e. a > 0.4854. 

We consider the transition between two spin network 
functions. Because of the chosen symmetry (|6|) within 
our lattice, this fixes all the spatial spins (and hence the 
initial and final spin networks to be equal). The only fluc¬ 
tuating spins are, therefore, the timelike ones, i.e. three 
per hypercuboid, six in total. Out of these three de¬ 
grees of freedom per hyper cube, there are two which are 
non-geometric. We disregard these in the following and 
observe that only two degrees of freedom remain in total 
for the two hyper cuboids. While one is obviously related 
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to the total 4-volume, the other corresponds to a (gener¬ 
alized) diffeomorphism direction connected to a time-like 
translation of the bulk Cauchy hypersurface (as in figure 
[2|). In other words, when computing the physical inner 
product between two spin networks, using the 2-complex 
described above, the sum over spins will, essentially, be 
the sum over two parameters, corresponding to the total 
4-volume of the universe, and a (quasi-)diffeomorphism 
degree of freedom. 

To shed light on the behavior of the amplitude with 
respect to this diffeomorphism, we consider the total am¬ 
plitude for the whole lattice, which is the product of two 
vertex amplitudes. The spatial spins are all set to some 
large spin j, while the time-like spin of the two hyper¬ 
cuboids are set to j(l + x) and jf(l — x), respectively. The 
amplitudes are normalized for x = 0. The findings are 
shown in figure [71 for different values of a. 



FIG. 7. Value of the product of two amplitudes I{a,x) de¬ 
pending on vertex translation x, for different values of the 
coupling constant a. 

Due to symmetry, it is clear that the derivative of the 
total amplitude I{a^x) = Av{j^j{l + x))Av{j,j{l — x)) 
is zero at x = 0. The second derivative depends linearly 
on Of, and obeys 


r'(«,0) = ^ - 12«- (29) 

This can be rephrased the following way: The contribu¬ 
tions from different points along the diffeomorphism orbit 
parametrized by x are different, and depend on a. There 
is a critical value ac = ~ 0.6067, and for a > ac^ 

the main contribution to the path integral comes from 
the point x = 0, i.e. the point where both hypercuboids 
are regular and of equal size. Conversely, for a < ac^ the 
main contribution to the path integral comes from very 
irregular configurations, i.e. the points x = ±1. Note 
that in the case of the path integral for more than two hy¬ 
percuboids, the situation will be even more pronounced, 
because the total amplitude factorizes over the vertices. 

One can therefore assert the following: the param¬ 
eter a distinguishes between two different regions, in 


which the path integral behaves quantitatively different. 
For a > ac, the main contribution to the path inte¬ 
gral, among all geometric configurations of equal volume, 
comes from the completely regular configuration, i.e. all 
spins, and therefore all hypercuboids, being equal. This 
can be interpreted as very regular geometry, approximat¬ 
ing 4-dimensional Euclidean space. 

On the other hand, for a < the main contribution 
to the path integral (again, among all geometric config¬ 
urations of equal 4-volume) come from configurations in 
which almost all 4-volume is concentrated in one time- 
slice, while all other time-slices contain nearly none of it. 
These configurations can be interpreted as very irregular. 

A comment is in order here: While this configuration 
can be regarded as higly irregular, it is nevertheless flat, 
and describes the same classical 4-d metric as the one in 
which all j’s are equal, i.e. the completely regular case. 
Although the j-values are different, these two states can 
be regarded as diffeomorphically equivalent, since they 
are results of a vertex translation symmetry as described 
in section Uni Geometrically, they correspond to different 
ways in which a 4-torus can be cut up into two hyper¬ 
cuboids - one in which the two hypercuboids have the 
same size, and one in which one contains much more vol¬ 
ume than the other. 

In other words, we have a clear example of a coupling 
constant a, which distinguishes between two regions in 
the phase diagram in which the path integral behaves 
qualitatively very different. Furthermore, on the criti¬ 
cal value a = ac, the amplitude function is nearly in¬ 
variant under change of x, which one can readily inter¬ 
pret as invariance under vertex translation symmetry, or 
diffeomorphism-invariance. This fosters the scenario for 
this point as second order phase transition: one can as¬ 
sume that correlation lengths become infinite (in terms 
of lattice distance) because of diffeomorphism invariance. 
The reason for this is that close and far separation on the 
lattice become diffeomorphically equivalent (see also dis¬ 
cussion in [illsl). Of course, this point deserves much 
more investigation, in order to make this correspondence 
more precise. 

One should note that the precise numerical value of 
ac should not be taken too seriously at this point: it 
rests on the simplifications we have employed, as well as 
taking all spatial spins to be equal. For the spatial spins 
different, one finds qualitatively similar results, however, 
with slightly changed values of ac. 


VII. PHYSICAL INNER PRODUCT AND 
TORSION 


A. The role of the physical norm 


The path integral can be used to compute the phys¬ 
ical inner product between states, and in particular to 
compute the physical norm of kinematical states. More 
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precisely, the path integral functions as a rigging map 

T] : ?^kin 3 I^kin ^ I>phys C T^phys (30) 

from (a dense subset of) the kinematical to (a dense sub¬ 
set of) the physical Hilbert space. Usually, Pkin is taken 
to be the linear span of the spin network functions. Dif¬ 
ferent spin network functions (which are all normalized 
in the kinematical inner product) will have images under 
the rigging map r] with different norm in Hphys- 

If a spin network state (j) on the (say, initial) boundary 
is given by spins je (and uniquely determined quantum 
cuboid intertwiners) for boundary edges e, then the phys¬ 
ical norm of such a state is given by 

( 77 ( 0 ) |?7(0))phys ^ n (31) 

jf V feedV 

where /e is the unique face whose intersection with the 
boundary is e. Note that in the boundary of the 
manifold has two connected components: the initial and 
the final spatial hypersurface. In an abuse of notation, (j) 
is considered to be a state on both. 


Ml M2 

FIG. 8. The separation of space-time M — 

This physical norm indicates how much the respective 
spin network contributes to the overall path integral, as 
we argue in the following. 

Consider a manifold A4 with boundary U Hj, being 
the product A4 = A4i#A42 of two manifolds in the sense 
of cobordisms (see figure [8]). If the Hilbert spaces are 
all finite-dimensional, then Z provides a functor in the 
cobordism category, satisfying 

Z(M) = Z(Mi) Z(M2). (32) 

For any observable O, its matrix element in boundary 
states = 77 ( 0 *’-^) therefore satisfies 

{¥\0\^f) = J2mv{<l^n)Uys{v{<Pn)\0\^%i.ys. 

(33) 

Note that in this last equation, the sum is being per¬ 
formed over an orthonomal basis {0n} in the kinemat¬ 
ical Hilbert space Hkin- By using the Cauchy-Schwarz 
inequality, each term in that sum is dominated by a con¬ 
stant times ||^(0n)||phys' other words, it can be ex¬ 
pected that terms with a very small physical norm do 
contribute very little to the matrix element. 


This argument can be generalized for arbitrary split¬ 
ting of the manifold: The two-complex F embedded in 
A4 is separated into two two-complexes Fi 7 ^F 2 . Then, 
each term in the state sum for Z(A4), i.e. each collection 
of spins (and intertwiners) on F, can be found as part 
of a term in the sum (|33|) . In other words, in the gen¬ 
eral state sum for Z(A4), one can approximately neglect 
the sum over all states which, for some decomposition 
(|32j) would result in a boundary spin network with small 
physical norm. 

Therefore, it is of interest for the evaluation of the 
state sum, to identify boundary states with small physical 
norm. By excluding the corresponding states from the 
state sum, one can drastically simplify the state sum, 
making only small approximation errors. In particular, 
the error can be estimated by the physical norm of the 
excluded state. 

It should be noted that the whole argument has to be 
slightly refined in the case of infinite-dimensional bound¬ 
ary Hilbert spaces. In that case, the state sum Z{M) 
does not exist as an operator on Hkin, but at most as a 
quadratic form on the dense subset Pkin- In particular, 
an equation like (|33|) will not hold in a straightforward 
sense. Still, since in this case the state sum will in many 
cases be defined as a limit of expressions which satisfy 
some sort of (|33|) . one can still expect the states with 
small (in particular those with vanishing) physical norm 
to be negligible. 

B. Torsion of boundary spin networks 

Motivated by this philosophy, we now compute the 
physical norm of boundary states for the special case of 
M ^ [0,1] X T^. Using the geometric interpretation used 
in section ini we will show that one can define boundary 
states which are non-geometric in the sense that they con¬ 
tain non-vanishing torsion. Furthermore, we will show 
that boundary states with large torsion have small phys¬ 
ical norm. 

In the continuum, the torsion tensor T is defined via 
— [ 65 , 60 ]^. Geometrically, T describes the 
failure of an infinitesimal square to close. After going 
around an infinitesimal parallelogram spanned by vectors 
and one gets translated by the vector 
(and rotated by the matrix 



FIG. 9. Torsion measures the failure of an infinitesimal par¬ 
allelogram to close. 

Taking literally the geometric interpretation of bound¬ 
ary geometries as provided by coherent polyhedra, this 
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suggests a way to identify torsion in the boundary spin 
networks we use in our analysis. 

The spatial hypersurfaces on the boundary are sepa¬ 
rated into cuboids, the areas of which are given in terms 
of the spins of the spin network function dual to it. This 
gives a geometric meaning to each cuboid, since its edge 
lengths are uniquely determined by its areas. These 
cuboids do not necessarily fit together without curving 
them, since for any square two of them are meeting, only 
the two respective areas agree, not necessarily the indi¬ 
vidual edge lengths. 

Consider a path going along a minimal square in the 
boundary spin network function (say, in the x^-plane), 
i.e. passing subsequently through the centres of four 
cuboids all meeting at one edge of the cubulation. The 
flat geometry suggests that each pair of neighboring sides 
are orthogonal to each other. So opposing sides are par¬ 
allel, but do not necessarily have the same length. The 
difference between the two lengths (say, in x direction) 
are precisely given by the T^^-component of the torsion 
tensor (see figure [10]). 



FIG. 10. Due to the non-geometric properties of twisted (spin¬ 
ning) geometries, the boundary geometry can contain torsion. 
In this case, T^y = -\- d 2 — — d^). 

With this in mind, we can construct spin networks 
which describe boundary geometry with torsion. For 
this, we choose a boundary of 2 x 2 x 2 cuboids. The 
spin network under consideration will have one (large) 
j distributed among the edges, apart from those on 
squares lying in the x^-plane. There are four of those, 
and they will each have one of the four possibilities of 
jf(l± A^) (see figure [TT]) . It can be readily seen that 
the components of the torsion tensor are T^y = 2Aa,, and 
T4 = 2A,. 


C. Results 

In our analysis, we have computed the physical norm 
m of the state '0 a^,a^ by using a hypercubic lattice of 

I 


2 X 2 X 2 X A, resulting in a norm 

\\i^Aa,,Ay\\‘N,a — ^ ( 34 ) 

timelikejf f v 

j(l + - Ay) 



FIG. 11. Boundary spin network with torsion: the four verti¬ 
cal lines (dual to the xy-squares in the boundary cubulation) 
carry spins j(l ± A^ ± A^), while all horizontal lines (dual to 
the xz and yz squares) carry the spin j. 


where the sum is being performed over all spins on time¬ 
like squares, while the spins over space-like squares are 
fixed by the lattice symmetry condition m- This means 
that the physical norm satisfies 

\\^A^,Ay\\N,a = (35) 


SO we can restrict ourselves to A = 1. In this case, there 
are six free spins in total, over which we integrate nu¬ 
merically, using Mathematica 10. We show the results, 
depending on the torsion of the boundary state, which is 
parametrised by A^,, A^, for three different values of a. 
The parameters Ax, Ay are varied from —0.4 to 0.4 in 
steps of 0.1, respectively. The results are normalized so 
that the state with A^, = Ay = 0 has physical norm one. 

The numerical calculation of the physical norm shows 
very similar features for all three values of a: The state 
with Ax = Ay = 0 have maximal physical norm (normal¬ 
ized to 1), while states with higher torsion have smaller 
physical norm. The norm of a kinematical state with 
torsion T is suppressed with a factor 

llV’TlIphys (36) 
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FIG. 12. Physical norm of states, depending on their torsion (color-coded). The plots are given for (from left to right) a = 4, 
a = 0.8, and a = 0.6. 


The rate of suppression Ca depends on a, and the con¬ 
nection between the two is depicted in figure [131 As one 
can see, the larger a, the more one can disregard torsion. 




0.8 1.0 


1.6 1.8 


FIG. 13. The inverse widths of the Gaussians in figure fT^ de¬ 
pending on the coupling constant a. Numerical investigations 
show that Ca = 17.220 — 10.45 is an excellent fit. 


It should be noted that our definition of torsion is a 
non-local property, in the sense that it requires more than 
one vertex to define it. It is closely related to the non- 
geometricity in the sense that torsion can only exist if 
there is non-geometricity in the vertices. However, the 
converse is not true: One can show that there are semi- 
classical configurations which are highly non-geometric, 
but have vanishing torsion. In that sense torsion-freedom 
is a weaker condition than geometricity. 


VIII. SUMMARY, CONCLUSION, AND 
OUTLOOK 

In this article, we have investigated properties of the 
quantum gravity path integral as given by the EPRL-FK 
spin foam model. We used the method of symmetry re¬ 
striction, which resulted in considering the sum only over 
spins and intertwiners which satisfy certain symmetry re¬ 
quirements. 

The discretization of space-time we considered con¬ 
sisted of a 4-dimensional hypercubic lattice, embedded 
in ^ [ 0 , 1 ] X R^. The intertwiners were restricted to 
the set of quantum cuboids^ i.e. coherent polyhedra with 
six faces, where the pairs of normals were negative of 
each others. 

Due to this choice, every edge was six-valent, and op¬ 
posite faces had to have equal spins. This enforced a large 
amount of symmetry throughout the lattice, restricting 
the configurations immensely. Despite the drastic simpli¬ 
fication, we gained several insights into the path integral, 
which should be useful in future analyses. 

The Barbero-Immirzi parameter 7 proved to be ir¬ 
relevant for the semiclassical analysis, which was a di¬ 
rect consequence of the fact that only flat conflgurations 
were considered. For the same reason, Newton’s constant 
could be set to 1 , and did not appear in the analysis. 

The only relevant coupling constant in our analysis was 
a, the power of the face amplitude. Several properties of 
the path integral crucially depend on it. 


Non-geometricity is suppressed in the path integral 

The states in spin foam models have a certain descrip¬ 
tion in terms of twisted (or spinning) geometries. It is 
well-known that these do not necessarily correspond to 
four-dimensional metrics. We identified the constraints 
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that states had to satisfy to describe a four-dimensional 
metric geometry, in terms of local constraints on the spins 
(i.e. areas). The constraints were equivalent to the exis¬ 
tence of an unambiguously defined 4-volume of the ver¬ 
tex. 

Deviation from this constraints were identified as non- 
geometricity Its physical role in the path integral cru¬ 
cially depends on the value of the parameter a. For 
a < ^ 0.49, out of all hypercuboids of a fixed 4- 

volume, the non-geometric ones contributed the most. 
On the other hand, for a > the geometric configu¬ 
rations dominate the non-geometric ones. In particular, 
for large enough values of a, non-geometricity is duely 
suppressed (see also the analysis in [53|). 

Furthermore, the effective action for the state sum in 
the large j limit contains a mass term m| for the non- 
geometricity which grows linearly with a. 

All of this analysis results from investigations around 
a regular hyper cuboid - for irregular hyper cuboids, as 
far as we can see, the qualitative features of the analy¬ 
sis remain the same, while e.g. the numerical value of 
changes slightly. As a general rule, it seems that 
non-geometricity is suppressed in the path integral, the 
stronger the larger a is. 

DifFeomorphism symmetry in the path integral 

As it turned out, despite the simplification of the path 
integral due to the restriction to quantum cuboids, the 
remaining state space is large enough so that an Abelian 
subgroup of the diffeomorphisms acts on it. In the ge¬ 
ometric large j-limit, this action can be interpreted as 
vertex-translation symmetry. It can be defined even 
on non-geometric configurations. Some of them, called 
semi-geometric, were such that the diffeomorphism ac¬ 
tion could be defined in a volume preserving way. The 
states could therefore be separated into three categories 

C semi- ^ geometric, 
states geometric ^ 

As it turned out, the value of the coupling constant a 
crucially influences the behaviour of the amplitude un¬ 
der the action of the diffeomorphisms. There is a critical 
value ac ~ 0.61, which separates two regions in phase 
space. The region a < ac is such that those configura¬ 
tions dominate the path integral, in which almost all of 
the 4-volume is concentrated in one time-slice, while the 
other time-slices does contain almost no 4-volume at all. 
On the other hand, for a > ac, the state with the major 
contribution to the path integral is the one in which the 
4-volume is separated equally among all vertices. 

A state (or a collection of states) which dominates the 
state sum can be interpreted as “vacuum”, in accordance 
with lattice gauge theory nomenclature. The reason is 
that it corresponds to ground states of the Hamiltonian, 
and the state sum describes fluctuations around it. We 
adopt this notion in the following discussion, and call 


the state (or the collection of them) which dominates the 
path integral as vacuum. 

The vacuum state in the two phases is quite different: 
while below the critical point a = o^c it is a superposition 
of space-times with one “fat time-slice” and many “thin 
time-slices”, above the critical point the vacuum is de¬ 
scribed by a regular lattice, i.e. where all spins through¬ 
out the hyper cubic lattice are equal. 

These two different states can be viewed as two differ¬ 
ent realizations of flat Euclidean space (i.e. a flat geom¬ 
etry on a torus). In one the vacuum describes the hy¬ 
per cubic lattice as very regularly embedded in the torus, 
while the other is a superposition of many different irreg¬ 
ular ways in which the hyper cubic lattice is embedded. 
Equivalently, the a > ac case can be seen as the com¬ 
pletely regular metric d(/)? (where ,... 04 are 

four angular coordinates on the 4-torus), while the highly 
irregular configurations dominating in a < ac can be re¬ 
garded as metrics of the form ds^ = 
some highly flutuating functions /i(0i),..., 04(04). Ob¬ 
viously, these two types of metrics are diffeomorphism 
equivalent. 

So the point a = ac separates two different vacua, 
both of which break diffeomorphism symmetry (i.e. the 
state-sum picks out different representatives of the dif¬ 
feomorphism orbit). 

There is furthermore some reason to assume that the 
critical point a = ac could constitute a phase transition 
of second order. The reason is that, at this point, ver¬ 
tex translation symmetry is nearly realized: to a certain 
numerical accuracy, all diffeomorhism-equivalent metrics 
describing the same geometry contribute the same value 
to the path integral. If one interprets vertex trans¬ 
lation symmetry as coming from diffeomorphisms, this 
would indicate that here one could find a diffeomorphism- 
invariant fixed point, even in the full path integral. 
Since diffeomorphism-invariance means that correlation 
lengths diverge (in terms of lattice distance), this could 
also indicate the existence of a continuum limit at this 
critical point. Obviously, this deserves a lot more inves¬ 
tigation, and is an exciting prospect for the search of the 
theory of quantum gravity. 


Torsion and the physical inner product 

In the last part of the article, we investigated the phys¬ 
ical inner product as defined by the symmetry-restricted 
state sum. In particular, we looked at the dependence 
of the physical norm of kinematical states, depending on 
the amount of torsion that was contained in its bound¬ 
ary geometry. We defined the torsion in terms of the 
naive geometric interpretation of the coherent polyhedra 
as cuboids in the large j-limit. It should be noted that, 
at this point, the connection to the quantum continuum 
torsion usually discussed in the literature mil is not 
completely understood at this point. This certainly war¬ 
rants further study. 
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What we found was that kinematical boundary states 
with high torsion have suppressed physical norm. The 
rate of suppression grows with increasing a. Since tor¬ 
sion, in the way that we have defined it, on the boundary 
implies non-geometricity in the bulk, this fits together 
with the observation that, for higher values of a, non- 
geometricity is suppressed in the path integral. However, 
the statement is stronger than that. The physical norm 
is not just the value of the amplitude for a specific state, 
but includes a sum over states. 

This has an important implication: To a good approx¬ 
imation, one can restrict the state sum to only those 


states, which are such that if the 2-complex is divided 
along any Sd hypersurface, the determined boundary ge¬ 
ometry on that hypersurface has no torsion. This is a 
non-local condition on states, given the naive definition 
of torsion we have used, in the sense that it is not a 
restriction of spins and intertwiners or vertices, but a 
property of the labels of bunches of neighboring vertices. 
Understanding this condition better, to find an effective 
approximation and restriction of the state sum, would be 
very interesting. It could also allow to get a much better 
and efficient way of computing state sums, at least for 
not too small values of a. 


^ 0.49 Q^c ~ 0.61 aoo = 0.75 

- 99 - • - 


non-geometric 
configurations dominate 


irregular configurations dominate 


geometric configurations dominate 


regular configurations dominate 


state sum finite 
for large lattices 


torsion suppressed in the physical Hilbert space 

-► 


FIG. 14. Different phases of the EPRL-FK model, depending on the coupling constant a. It should be noted that the precise 
numerical values should not be taken too literally, due to the simplifications and choices we have made in our analysis. However, 
we do expect that a qualitatively similar phase diagram could be found also in the full path integral. 


Outlook 

Although some interesting statements about the spin 
foam state sum could be derived in this article, it re¬ 
mains to be shown that these statements survive if more 
states are included than just the ones containing quan¬ 
tum cuboids. In particular, it is paramount to include 
states which allow for the summation over curvature de¬ 
grees of freedom, to arrive at a more realistic approxima¬ 
tion to the whole state sum. 

Furthermore, it should be noted that we used the un¬ 
renormalized EPRL-FK amplitude. In working on a fixed 


lattice, this disregards radiative corrections from finer 
lattices. 0 What one should actually do is use the renor¬ 
malized amplitude instead, which serves as an effective 
action on the fixed lattice. This should in particular con¬ 
tain more coupling constants, such as e.g. the cosmolog¬ 
ical constant or prefactors in front of higher curvature 
and even non-local terms. 

At this moment, however, this amplitude does not ex- 


In the GFT picture, this would correspond to sum over more 
than one triangulation. 
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ist. In fact, the situation is such that very little is known 
about the renormalization group flow of the EPRL-FK 
model. Performing actual calculations in this direction 
with the full state sum might prove very challenging, both 
analytically and numerically. This is why, to tackle this 
problem, the methods shown in this article, could prove 


useful to that end. 
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